Adiabatic dynamics in open quantum critical many-body systems 
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The purpose of this work is to understand the effect of an external environment on the adiabatic 
dynamics of a quantum critical system. By means of scaling arguments we derive a general expres- 
sion for the density of excitations produced in the quench as a function of its velocity and of the 
temperature of the bath. We corroborate the scaling analysis by explicitly solving the case of a 
one-dimensional quantum Ising model coupled to an Ohmic bath. 



The understanding of the non-equilibrium dynamics 
of strongly correlated quantum systems is one of the 
most challenging problems of modern condensed mat- 
ter physics. Interest on this subject has been revived 
recently by unprecedented experimental breakthroughs 
in the context of cold atomic gases (see e.g Non- 
equilibrium conditions in cold gases can be realized con- 
troUably in various ways, e.g. by a proper choice of the 
initial state, or by changing the Hamiltonian in time. 
In view of the variety of situations that can be studied 
in this context, it is of paramount importance to find 
paradigmatic situations that allow the general features 
of the non-equilibrium dynamics of many-body systems 
to be understood. One such paradigm is obtained when 
the parameters of a quantum system close to a quan- 
tum phase transition are varied in time in such a way as 
to traverse the quantum critical point. Because of the 
vanishing of the gap A at criticality, a finite density of 
defects is generated, no matter how slow is the quench, 
as first shown in Ref. [2, 3]. The density of defects is a 
universal scaling function of the quench velocity v, as in 
the Kibble-Zurek (KZ) mechanism Q originally derived 
for classical continuous phase transition. In addition to 
its intrinsic interest, this problem is relevant to adiabatic 
quantum computation |5j| and quantum annealing 

The intense theoretical activity following Refs. 0, S| 
has clarified several important issues (see |7| and refer- 
ences therein) on the adiabatic, phase coherent dynamics 
of closed many-body critical systems. Closed systems, 
however, are only idealizations: any quantum system is 
weakly coupled to an environment inducing relaxation 
and dephasing. This observation motivated a series of 
recent studies, in particular on the effect of classical [9] 
and quantum [lo| noise acting uniformly on a quantum 
Ising chain, and on the effect of local noise on a disordered 
Ising chain with up to 20 spins [ll|. However, the most 
natural and important question remained answered: to 
what extent is it possible to describe universally the pro- 
duction of defects in an adiabatic quench in the presence 
of dissipation and dephasing ? In this Letter, we airswer 



Figure 1: A sketch of the crossover diagram at finite tempera- 
tures due to the presence of the quantum critical point, where 
the gap A — » 0. The system traverses the quantum critical 
region in a time tgc- 



this question by showing with a general scaling analysis 
that even in the presence of an external environment the 
adiabatic dynamics of open critical systems is governed 
by universal scaling laws (with modified exponents). To 
support this statement, we solve by means of quantum 
kinetic equations the adiabatic dynamics of a quantum 
Ising chain coupled to a local external environment, lead- 
ing to the relaxation of all quasiparticle modes. 

Scaling analysis of defect formation. At equilibrium, 
all the characteristic features of a quantum phase transi- 
tion leave an imprint at low temperatures, leading, close 
to the quantum critical point, to a crossover at temper- 
atures T '--^ A. For T <C A, the low frequency dynamics 
can be described semiclassically in terms of quasi-particle 
excitations, while in the quantum critical region (red re- 
gion in Fig[T|^T A and quasi-particle excitations no 
longer exist [8| . The purpose of the analysis below is to 
discuss the universal fingerprints left on the nonequilib- 
rium dynamics by the quantum phase transition in the 
presence of a bath. For this sake, we now start by pre- 
senting a scaling analysis describing the influence of a 
bath at temperature T on the production of defects. 

We start by considering a linear quench of the control 
parameter h from an initial value hi to a final one hf 
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across the critical point he at a speed v. The system is 
initially in equilibrium with the bath, whose temperature 
T is kept constant during the quench. This is sketched in 
FigU The passage through the quantum critical region 
(T :s> A) leads to substantial heating effects because of 
the relatively high bath temperature with respect to the 
characteristic energy scale A for the system excitations. 
We shall see that, as long as this is the most important 
time interval for the bath-systems interaction, only the 
low energy details of the system spectrum matter, and 
universality holds. 

The universality of the defect production in the pres- 
ence of a bath follows from two conditions, discussed be- 
low. First of all, the density of excitations £ can be 
written as 

£ — £kZ + £inc (1) 

where £kz is the coherent contribution, present also in a 
closed system, and Sine is the incoherent one due to the 
presence of the bath. For a closed system, the density of 
excitations was found to scale as Ekz oc v'^^^'^^"^^^ 
where d is the spatial dimension of the system and z 
are the correlation length and the dynamical critical ex- 
ponents respectively. Moreover we assume that the bath 
does not influence the system in the semiclassical regions 
(T ^ A). Hence, we consider thermal excitations pre- 
dominantly created inside the quantum critical region, 
where T ^ A, at a rate c>c aT^, a being the 

system-bath coupling constant. The dynamics for the 
population of the (excited) mode k is thus governed by 
the rate equation ^Pk — —r^^iPk — Pl*^ {he)), where 
Pf^ihc) is the critical thermal equilibrium distribution 
inside the quantum critical region. Integrating this rate 
equation we obtain Pk{hf) - (1 - e-^"'*«^)P^'' (/i^), 
where tqc — 2T^/'^^t;"^ is the time spent going through 
the quantum critical region (see FiglT]), whose bound- 
aries are given by T ~ A ^ \h — hc^^^ . Integrating over 
all fc-modes we obtain: 

£^nc oc (l - e-^) j dE Ei-'Pl'^ (K) , (2) 

where we used the scaling of the critical energy E oc k^. 
Finally, since the thermal distribution P^'' is a function 
of E/T, changing variable to E/T leads to 

£„,c (X av~^T'>+^ , (3) 

vahd in the limit T-^^"^^ <^ vt. The scaling of Sine 
is directly related to the time spent inside the quantum 
critical region. The crossover from the coherent to the in- 
coherent defect production is reached when Sine — £kz, 
giving 

Vcross OC a -(-"wi t(^+ ) . (4) 

The different scaling with v of £kz and £i„c implies that 
for fast quenches, v > v cross, the KZ contribution domi- 



nates, while for slower quenches, v < Vcross, the incoher- 
ent contribution due to the thermalization induced by 
the bath is the most important. 

Eqs. (IT|),([2]),(|11) represent the generalization of the scal- 
ing laws given in Refs.[l, Q to the case of open quan- 
tum critical systems. They are in principle amenable of 
an experimental veriflcation and constitute the key re- 
sult of this work. In the case of a quantum Ising chain, 
which we discuss in detail in the second part of the Let- 
ter, the previous expressions specialize as follows. The 
time spent within the quantum critical region scales as 
tqc — 2Tv~^, since = z = 1. We consider an Ohmic 
bath that acts as a random external magnetic field on 
each site (see Eq. ([S])). We find that the relaxation time 
in the quantum critical region scales as r^^ ~ aT^ (see 
Fig. [5]), i.e., 9 — 2, in agreement with a Fermi golden 
rule argument. It follows that the contribution to the 
defect production induced by the bath and the crossover 
velocity scale as Sine oc av~^T'^ and Vcross oc a^/'^T^/^, 
respectively (see Fig. [3]). 

Quantum Ising Model and kinetic equations. To sup- 
port Eqs. ^ and ^ we now study the physics of a quan- 
tum quench for a quantum Ising model coupled locally 
to a set of Ohmic baths. The locality of the system-bath 
coupling causes the breaking of the translational symme- 
try of the closed system, hence allowing the discussion of 
the quench dynamics in the presence of relaxation of all 
elementary excitations. Notice that no qualitative fea- 
tures are expected to emerge in the case where the baths 
are correlated over a finite distance because (in the scal- 
ing limit) the correlation length is the largest lengthscale 
in the problem, and details of bath correlations over mi- 
croscopic distances should not matter. It is also impor- 
tant to observe that long time correlations induced by 
the bath may change the universality class of the tran- 
sition 15]. Despite the intrinsic interest of this issue, 
which found recently application in the context of the 
physics of cold atoms J^i, we will not consider it here. 
Therefore, we will further assume that the bosons have a 
non-zero inverse lifetime 7 ^ T which provides a natural 
cutoff-time for the bath correlation functions. 

The model we consider is then defined by the Hamil- 
tonian 

7 ^ 

^ = -2E{^J^^+i + [^W + ^^]'^|}+^^ • (5) 

j 

It consists of a chain of N spins {a^ and cr^ are Pauli 
matrices) with an Ising interaction and subject to a 
transverse magnetic field h{t). The bath couples to cr^, 

with Xj = Z]/3 ^/3(^/3,i + bpd), where b'^^ jibp.j) are the 
creation (annihilation) operators for the bosonic bath 
modes coupled to the j-th spin. The bath Hamilto- 
nian is Hb /B^^l^^j^j^l^j- The system-bath coupling 

is chosen to have Ohmic spectral densities A^(5(w — 
LUf}) = 2auj eKp^—Lo/uc), where lUc is a high-energy cut- 
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off 13|, [iJI . In the case of no coupling to the bath {a = 0) 
the system has a QPT at /ic = 1, and for h < he a 
spontaneous magnetization along x appears. The gap 
A = \h — h(.\ induces at finite temperature a V-shaped 
crossover phase diagram @ , as sketched in Fig. [TJ 

We now analyze the problem by deriving a quantum 
kinetic equation which allows us to calculate the density 
of defects produced after the quench. This procedure al- 
lows us to describe the effect of the environment also in 
regimes which are beyond the applicability of the scaling 
laws deduced above, where universality is not expected 
to hold. In order to describe the dynamics of ^ it is 
first convenient to map the spins onto spinless fermions 
by means of a Jordan- Wigner transformation. In mo- 
mentum space, the Hamiltonian ([5]) reads: 



H 



k>0 



1 



— y^^lf^^^k+,X, + HB, (6) 



where Hk = — [cos {k)+h{t)]Tz+sm {k)Ty, f are Pauli ma- 
trices in the Nambu space defined by the two-component 
fermion '1'^ = ( c\ c_fc ). When Xq — 0, the subspaces 
of the different k modes are decoupled, the Hamiltonian 
is quadratic and can be diagonalized by a Bogoliubov 
transformation yielding H = J2k>o ^^ivlvk - V-kV-k) 
with Ak = Vl + 2/i cos k + h'^. The interaction with the 
bath causes a mixing of modes with different momenta. 
We determine the density of excitations by deriving a 
quantum kinetic equation for the Green's function using 
the Keldysh technique [l^. This kinetic equation will 
be expressed in terms of the fermion lesser Green's func- 
tion , which is a 2 X 2 matrix in Nambu space with 
components -i[G<{t,t)],^j = Using a 

self-consistent Born approximation for the bath-mediated 
scattering of k-modes, valid for weak bath-system cou- 
pling a <C 1, and the Markov approximation (justified by 
the assumption of a cutoff time for the bosonic modes), 
we get: 



dtG< + i nk,G 



(7) 



lY,r^{i + ^G<)D,kG< + G< olil + zG<) + H.c. 



N 



where we neglect irrelevant Lamb shifts. Here Dqk = 
i j^dsg>{s)U%t ~ s)r-Ukit,t - s), where g>{t) = 
— i {Xq{t)Xq{Q)), and Uk{to,t) is the evolution operator 
satisfying idtUk — 'Hk{t)Uk- Parameterizing —iG'^ = 
1/2 [i - (1 - 2Pk)fz + CkT+ + Glf-] after the Bogoli- 
ubov transformation which diagonalizes 'Hk^ one finds 
that Pk = {vivk) and Gk = (ij-kilk)- Therefore, the 
density of defects produced can be expressed as 



£ = 



2N 



ETr[(i 



1 



fc>0 



(8) 
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Figure 2: r ^ /T'^ as a function of A/T (negative values cor- 
respond to h < he). Data are relative to different distances 
S = y/T^ + {h-h^y from the QPT. Close to the critical 
point [5 0) all curves collapse into a unique scaling func- 
tion. Dashed line is a fit a(l -I- bA/T) exp {-A/T} with 
a ~ 9.4(6) and 6 ~ 0.9(3). Inset shows the data collapse 
of T~^/A^ as a function of A/T for different values of h: 



oc exp{-A/r} for T < A. 



By solving the kinetic equation ([7]) we are able to calcu- 
late £ in Eq. ([8]), and thus to analyze the Kibble-Zurek 
mechanism in an Ising chain coupled to a bath. The 
results of this analysis are presented below. 

Kinetic equation and the scaling regime. As dis- 
cussed in the first part of this Letter, in order to de- 
termine the scaling at finite temperature we need to 
know the relaxation time r for the excitations. To 
this end, it is sufficient to consider only the dynamics 
of the populations Pk in the kinetic equation and ne- 
glect the off-diagonal components. By linearizing the 
kinetic equation around the equilibrium (Fermi) distri- 
bution, one gets dtSP = —TZ{h, T)SP, where SP = 



Pk, 



Tjth 
^kt ' 



,Pk 



N/2 



P, 



th 



The characteristic re- 



laxation times of the system are the inverse eigenvalues 
of TZ. The leading asymptotics at long times for all pop- 
ulations Pk is governed by the inverse of the smallest 
eigenvalue Ai of 7?., t = Aj~^. By numerical inspection, 
we find that r'^ - aT^fi^) e'^'^ (see Fig [2]). 

The scaling obtained for the relaxation time, together 
with an explicit integration of Eq. ([2]) for the Ising model, 
lead to 



log 2 
27r 



T 1 



-2T/{tv) 



(9) 



fc>0 



which, in the limit 2T/tv <^ 1, confirms the scaling 
result given in Eq. ([3]). In Fig. [3] the solution of the 
kinetic equation is compared with the ansatz given in 
Eqs. H]) and ([9]). The agreement is excellent, confirming 
our scaling approach. The crossover value of v which sig- 
nals the transition from the coherent- to the incoherent- 
dominated defect production obeys the power-law scaling 
given by Eq. At lower quench rates the full expres- 
sion, Eq. ([5]), is needed for an accurate comparison with 



4 




Figure 3: Left: density of excitations E versus the quench ve- 
locity V. Circles are obtained by integrating the kinetic equa- 
tion (KE); solid and dashed lines are Eq. ([1} with fi„c given 
by (|9]) and its linearized expression, respectively. KZ scal- 
ing £kz oc \pv i2j] is plotted for comparison as dotted lines. 
Plots refer to iv'= 400 and a = 0.01 and T = 0.15, 0.1, 0.07 
(from top to bottom). The final values of magnetic field are 
hf — 0.8 (outside the quantum critical region for the val- 
ues of T considered) and the critical point hf — 1. For 
hf = 1 we considered tQc = v~^T since only half of the 
quantum critical region is crossed, a) Data collapse of £i„c 
obtained from the kinetic equations by considering only in- 
coherent thermal transitions (rhs of ((Zl); data refer to both 
hf = 0.8, 1 and lO""* < v < 10"^. The fit verifies the scal- 
ing ((3| Sine OC v~^T'^. In b) the scaling of Vcross IS obtained 
equating £i„c from kinetic equations to £kz for hf = 0.8, 1; 
the fit confirms the scaling predicted by Q, Vcross oc T^/'^. 




Figure 4: Same as left panel of Fig. O but for hf = 0.6; 
curves are relative to T = 0.15 only. Stars represent the 
density of excitations produced incoherently by the bath and 
are obtained by integrating the kinetic equations and ignoring 

the unitary evolution term i \^k, G^j , responsible for the 

coherent excitation process. The non-monotonous behavior 
is due to the non-critical relaxation mechanism, that is more 
relevant at low v. In the regime of slow v scaling no longer 
holds and only thermal excitations contribute to £. 



the solution of the kinetic equation. 

On lowering the final value hf oi the field, the agree- 
ment with the scaling ansatz becomes worse at low 
quench rates v (see Fig. HJ. This is due to a non-critical 
relaxation mechanism that depends strongly on the de- 
tails of the energy spectrum for hf < h < he- Once 
the system leaves the quantum critical region, entering 
the semiclassical region where T ^ A, the bath starts 



to relax-out the excitations previously created; if time 
spent in the semiclassical region is long enough all exci- 
tations disappear. Outside of the weak coupling regime, 
when the relaxation time and quench time are compara- 
ble, scaling is expected to be recovered in thermodynamic 
quantities, such as total energy, heat, entropy and work 
done on the system ITj . 
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